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In this paper, the authors develop a dynamic model of a commercial micro-combined heat and power 
(mCHP) unit and analyse its dynamic behaviour when the engine is running at different mass flow in¬ 
puts. The simulation predicts with a low root mean bias error (RMSE) the most important outputs from 
the cogeneration unit during the starting, steady-state and stopping periods. Furthermore, the presented 
transient model reproduces the behaviour of the cogeneration unit when the fuel and air mass flows are 
changing. The obtained results are discussed, and the different possibilities for the variation of the 
thermal to power ratio are analysed. These combinations include the variation of the flow distribution 
inside the machine and the position of the exhaust heat exchanger. The power to thermal ratio can be 
modified between 0.15 and 0.26 for these combinations. The performance of the engine and the variation 
of the heat source temperature are also analysed theoretically. The simulation results conclude that an 
important saving could be obtained when the electrical to thermal ratio (ETTR) is tracked for the power 
or thermal demands from a dwelling. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Stirling engines (SEs) will represent an important micro¬ 
cogeneration device in the next few years because classical en¬ 
ergy resources are progressively increasing in cost. The external 
combustion of SEs allows the utilisation of different heat sources, 
which makes SEs more independent of energy market fluctuations. 

Moreover the utilisation of mCHP units is expanding in the 
domestic sector because different technologies are being developed 


* Corresponding author. Tel.: +34 986 818 624. 

E-mail addresses: acacabelos@uvigo.es, antoncare@hotmail.com (A. Cacabelos), 
peguia@uvigo.es (P. Eguia), jmiguez@uvigo.es (J.L. Mlguez), guillermo.rey@cud. 
uvigo.es (G. Rey), elena.arce@uvigo.es (M.E. Arce). 

http://dx.doi.Org/10.1016/j.applthermaleng.2014.07.078 

1359-4311/© 2014 Elsevier Ltd. All rights reserved. 


[1 , and the emerging heat flow market represents an important 
role in district heating systems. Among others, the Stirling engine 
(SE), internal combustion engine (ICE), organic Rankine cycle (ORC) 
and fuel cell (FC) technologies have great potential for development 
[1 . The electric power output of these domestic devices is normally 
in the range of 0-2 kWe, although the electrical to thermal ratio 
(ETTR) varies between 0.08 and 0.59 depending on the technology. 
In SEs, this variation is still smaller. Normally, this ratio is reduced 
to the interval 0-0.3 in the mCHP unit with Stirling technology. 
Although cogeneration capabilities are the most common, some ICE 
devices can also supply the dwelling with cooling during the 
summer time using tri-generation capabilities 2,3 . 

There are different studies that simulate the operation of a SE 
cogeneration unit. Some studies may discretise the performance 
time into different running modes [4,5] to take into account the 
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transient behaviour. Other studies may integrate the cogeneration 
unit into a building simulation environment to analyse the coupling 
with a dwelling during a meteorological year [6 . Other works may 
focus on this coupling and may analyse different SE devices of the 
building sector for the seasonal performance and the desired 
characteristics from the cogeneration unit [7,8 . 

The upper limit of the Stirling engine efficiency is bounded by 
Carnot [9] and cannot be exceeded. Depending on the charac¬ 
teristics from the working fluid (principally conductivity, specific 
heat, viscosity and density) within the closed thermodynamic 
cycle and the working pressures and temperatures, the electrical 
to thermal ratio can be strongly modified, and power efficiencies 
between 10 and 60% could be obtained 10]. Normally, this ratio is 
constructively designed and depends principally on swept and 
dead volumes, on regenerator characteristics and on the selected 
working fluid [6]. However, these variables are difficult to change 
during the running mode. Some SEs can externally modify the 
stroke, the phase angle or the dead volume from the piston, but 
these capacities require an extra mechanical complexity, which 
for many instances is not justified. Moreover, there are some 
types of SEs where these modifications are impossible. Although 
the phase angle variation is one of the best methods for power 
control, the double-acting engine in this study does not allow this 
possibility 11 . 

The power control of the SEs can also be performed by varying 
the operational characteristics, such as pressure, temperature, 
speed or electric load. Occasionally, the type of fuel used or the 
engine configuration restricts the power control method 1 . Other 
authors analyse the behaviour of the SEs through the variation in 
the characteristics of the working fluid. This working fluid is 
simulated as a lumped-mass transient model. The shaft power 
output is investigated when different parameters from the SE are 
changed [12]. This model reveals that as the operating speed in¬ 
creases, the difference in temperature between the source and the 
sink decreases. Consequently, the shaft power output also declines, 
but a maximum is reached at a critical speed. 

A lump transient model is also a valid method that is often used 
to simulate the cold start modelling and transient response of 
thermal engines [13 . 

Stirling engines were also dynamically simulated using dynamic 
models based on experimental data [14] or by determining the 
parameters that the IEA/ECBS Annex 42 specifies 15], although 
these methods are sometimes cumbersome and difficult to 
undergo. 

Recent studies analyse the feeding of a family house by means of 
external and internal combustion engines (ICEs), and the results are 
compared with traditional power production systems [16]. The 
three selected micro-cogeneration devices were operated in heat 
manage mode to control the temperature of a thermal storage tank, 
and only the ETTR that was established by the manufacturers was 
used. The modification of this ratio becomes especially interesting 
in isolated applications, where a net connection or district heating 
is not available or when the excess energy is under-paid. 

Comparing ICEs with external combustion engines, which are 
the most common technologies in mCHP units, the ICEs present 
higher mechanical losses because of more moving parts and higher 
noise and vibrations due to the explosions and high operating 
speeds. ICEs were also calibrated using the IEA Annex 42 17] but as 
previously mentioned, this method is time-consuming and requires 
the estimation of more than one hundred coefficients, which are 
derived from experimental data. 

Most of the SE cogeneration units operate at full capacity and try 
to work at the most efficient working point. While the temperature 
from the heat source is kept as high as possible, the mean cycle 
pressure is also kept high and consequently the power output 


because they are directly proportional. Normally, there is a 
maximum power at a determined engine speed. This maximum 
will depend on the working pressure, source temperature, the 
working fluid or the regenerator wire mesh 18]. This work reveals 
that a lineal relation also exists between the engine speed and the 
obtained torque in the shaft. 

Thermal losses can have an enormous influence on the results. 
Although some analyses do not consider these losses, most of them 
include irreversibilities and thermal losses into the cooler, the 
heater and the regenerator 19,20]. Numerous works analyse the 
behaviour of different SE configurations based on different ther¬ 
modynamic assumptions [21,22 . Others analyses try to validate the 
theoretical models, which simulate the thermodynamic behaviour 
of an experimental SE unit [18,23]. 

Although numerous works have centred around the Whisper- 
gen unit [5,10], no analysis has been oriented to create easy use and 
to calibrate the model's integration into a simulation tool envi¬ 
ronment with a transient and a variation of the ETTR capacities. 

Therefore, the present paper describes and simulates a com¬ 
mercial SE design with a double-acting alpha configuration. A 
previous lump model from the authors [24 was improved because 
when this model was applied in transient test varying the mass 
airflow going into the combustion chamber and consequently the 
fuel consumption and source temperature, it was discover that this 
model did not perfectly follow the transient behaviour from the 
machine. Furthermore, the analysis that was undergone in this 
article includes a transient study of the cogeneration device and 
different possibilities of varying the ETTR during the running mode. 
Consequently, this ratio could be slightly modified attending to 
household demand or user's interests. That modulation does not 
include variations in the internal parts of the engine such as the 
material or porosity of the regenerator or modifications in SE's 
geometry. Cycle mathematical analyses were not considered and 
only variations in mass flows or source temperatures including 
different heat exchanger arrangement have been evaluated. 

Furthermore, this transient model could be applied in different 
scenarios especially in isolated applications such as the presented 
in Refs. 25] or 26 . The model is also interesting for grid connected 
or heat flow market applications where the variation of the ETTR 
could be applied depending on market prices or user's demands at 
any time. 

In the first sections of this article, the main characteristics of the 
selected Stirling engine are explained focussing on the power 
output at different source temperatures. Subsequently, the 
improved transient model based on two lumps together with the 
procedure to integrate it in the simulation environment TRNSYS are 
explained. An analysis about the main exchanger energy inside the 
cogeneration system is made then and the new proposed model is 
presented. Finally, the experimental and simulated results are 
discussed and the energy delivered during the shutdown process is 
analysed. 

2. Materials and methodology 

Firstly, different tests were carried out for a test bench princi¬ 
pally consisting of a commercial Stirling device from the manu¬ 
facturer Whispergen together with an electric circuit and a thermal 
circuit to simulate the hypothetical electric and thermal loads from 
a dwelling. 

These tests were performed in transient mode and were 
running on diesel and with different air mass flow inputs. The 
procedure was similar for all of the one-hour length tests. These 
data were used to calibrate the lumped mass model developed in 
the previous article by the authors in Ref. [24] and was modified to 
better fit the transient behaviour. In this way, the new dynamic 
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system presents a better transient response when the airflow and 
consequently, the fuel input entering into the burner is varied. 

Once the model is validated, the cogeneration unit is analysed, 
and the possibilities for regulating the electrical to thermal ratio 
(EUR) are investigated. The easiest way to achieve this regulation 
using the same working fluid and Stirling morphology is to adjust 
the hot and cold source temperatures, which directly modifies the 
maximum efficiency from Carnot's cycle [9]. The sink temperature 
is delimited by the type of heating system adopted, but the hot 
source temperature could be slightly changed for the loads or the 
consumer's interests. 

3. Problem definition 

The need to research cogeneration devices efficiency in its final 
location caused the emergence of transient models. Normally, these 
transient models require the calculation of numerous constants 
which involves time dedication. The present work proposes an 
easy-to-calibrate model which is perfectly integrated in a simula¬ 
tion environment (i.e. TRNSYS software) for subsequent simulation. 
Moreover the possibilities of variation of the ETTR from the 
cogeneration unit are investigated. Thus, one proposed scenario 
analyses position variation of the exhaust gases heat recovery 
exchanger. This modification incorporates important changes in the 
EUR which may be used to increase the efficient of the system in 
particular cases. 

Constant EUR is normally used in research because units are 
optimised to run within the maximum efficiency. However the 
cogeneration unit under study was developed to run in an auton¬ 
omous way; so the ETTR variation capacity becomes in a very 
interesting feature. 

4. Characteristics of the SE 

The selected SE is a commercial cogeneration device from the 
manufacturer Whisper Tech Limited. The model is the PPS16-24MD 
that corresponds with a double-acting alpha-type model, which 
uses a wobble yoke mechanism to transform the reciprocating 
movement of the pistons in a rotational one to move an alternator. 
This model burns diesel, but it can be easily transformed to run 
with a different fuel. Some projects made this engine work with 
bio-oil or ethanol 10], and the manufacturer actually commer¬ 
cialised a similar device that was fed by natural gas. 

The engine presents four cylinders with pressurised nitrogen at 
a maximum of 28 bars. Table 1 shows the principal SE parameters 
taken from Ref. 27 . The dead volumes are approximately an order 
of magnitude less than the swept volumes and although the pres¬ 
ence of dead volumes reduces the power output from the system, 
the efficiencies are not always reduced 11 . 

Fig. 1 represents the distribution of the different volumes and 
the heat exchangers inside the SE. The wobble yoke mechanism 
together with the alternator is also represented. The double-acting 
alpha-type model is a special arrangement where the lower part of 
each cylinder acts as the displacer of the adjacent piston as shown 
in Fig. 1. 

The experimental installation of the test bench, the data 
acquisition system together with sensors' distribution and time 
step selection were explained in a previous article from the authors 

of [26]. 

5. Analysis and simulation of the power output 

The Schmidt analysis is a common theoretical method where 
the different configurations of SEs were highly analysed. It was 
performance the Schmidt analysis into the Whispergen unit to 


Table 1 

Variables and Stirling engine parameters. 


Description 

Variables 

Value 

Units 

Specific heat of flue gases 

C P- gas 

1.107 

[kj/kg K] 

Specific heat of water 

Cp_ water 

4.18 

[kj/kg I<] 

Thermal capacitance of block 

MCblock 

13.625 

[kJ/K] 

Thermal capacitance of burner 

MCburner 

54.75 

[1<J/K] 

Air mass flow rate 

™air 


[kg/s] 

Diesel mass flow rate 

771 f U ei 


[kg/s] 

Flue gases mass flow rate 

ffigas 


[kg/s] 

Heating system water mass 

f^water 


[kg/s] 

flow rate 

Lower heating value 

LHV 

42,800 

[kj/kg] 

Energy produced with gas 

Qcomb 


[kj] 

oil combustion 

Energy absorbed by the water 

Qwater 


[kj] 

from the lumped mass 

Root mean square error 

RMSE 


[-] 

Ambient temperature 

Tamil) 


[K] 

Experimental temperature 

Texp 


[I<] 

Exhaust gases temperature 

Texh 


[K] 

Flue gases temperature 

Tgas 


[K] 

Temperature of lumped block mass 

Tblock 


[I<] 

Temperature of lumped burner mass 

Tburner 


[I<] 

Water temperature at heat 

Tout 


[I<] 

exchanger outlet 

Water temperature returning 

Tret 


[K] 

from heating system 

Simulated temperature 

Tsim 


[I<] 

Overall heat transfer coefficient 

(U-A)iumps 


[W/I<] 

from burner to water block 

Overall heat transfer coefficient 

(U-A)he 


[W/K] 

for gases-water heat exchanger 

Density of lumped mass 

P 


[kg/m 3 ] 

Bore 

B 

4.200E-02 

[m] 

Stroke 

S 

2.000E-02 

[m] 

Piston area 

A P 

1.385E-03 

[m 2 ] 

Swept volume expansion piston 

^SE 

3.414E-05 

[m 3 ] 

Swept volume compression piston 

^sc 

2.813E-05 

[m 3 ] 

Dead volume of expansion space 

Yde 

3.592E-06 

[m 3 ] 

Dead volume of compression space 

Vdc 

6.751 E-07 

[m 3 ] 

Heater volume 

V h 

9.180E-06 

[m 3 ] 

Cooler volume 

V k 

8.269E-06 

[m 3 ] 

Regenerator volume 

Vr 

1.745E-05 

[m 3 ] 

Total mass of working gas 

M 

3.750E-04 

[kg] 

Thermal resistance from 

Rh 


[m 2 K/W] 

the hot side 

Heat flow going into the 

Qin 


[W/m 2 ] 

expansion space 

Convection coefficient in the 

h H 


[W/m 2 K] 

expansion space 

Thermal resistance from 

R c 


[m 2 K/W] 

the cold side 

Heat flow going out the 

Qout 


[W/m 2 ] 

compression space 

Convection coefficient in the 

he 


[W/m 2 K] 

compression space 

Expansion space gas temperature 

Th 

9.750E+02 

[K] 

Compression space gas temperature 

Tc 

3.600E+02 

[K] 

Regenerator space gas temperature 

Tr 

6.675E+02 

[I<] 

Phase angle 

A 

9.000E+01 

[°] 

Temperature ratio 

T 

3.692E-01 

[-] 

Swept volume ratio 

V 

8.239E-01 

[-] 

Dead volume ratio 

X 

1.880E-01 

[-] 

Engine speed 

N 

1.500E+03 

[rpm] 

Engine speed 

N 

2.500E+01 

[Hz] 

Indicated expansion energy 

w E 

3.651 E+01 

[J] 

Indicated compression energy 

W c 

1.348E+01 

[J] 

Indicated energy 

Wi 

2.303E+01 

[J] 

Number of cylinders 

z 

4.000E+00 

[-] 

Specific heat at constant 

Cpn 

1.039E+03 

[J/kg I<] 

pressure [Nitrogen 300 K] 

Specific heat at constant 

r 

*~vn 

7.430E+02 

[J/kg I<] 

volume [Nitrogen 300 K] 

Gas constant 

R 

2.968E+02 

[J/kg K] 

Gas constant 

R 

8.314E+00 

[J/K mol] 

Vde/^se 

^DE 

1.052E-01 

[-] 
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Table 1 ( continued ) 


Description 

Variables 

Value 

Units 

Vdc/^SE 

XdC 

1.977E-02 

H 

Vr/Vse 

Xr 

5.110E-01 

H 

Molar mass N 2 

Mn 2 

2.801 E-02 

[kg/mol] 

Thermal conductivity N 2 

A 

2.598E-02 

[W/K m] 

Number of poles 

P 

12 

H 


better understand its behaviour when expansion temperature is 
varied. 

This mathematical model does not take into account the irre¬ 
versibility introduced by friction, throttling or pressure losses 
because it is a simple method where the working gas is considered 
an ideal gas. Although there are different improvements to this 
theory as introduced by Gustav, Finkelstein, Senft [1] or Petresku 
[28] that attempt to include the internal losses generated into the 
engine, our dynamic lump model tries to simulate the global 
behaviour of the SE in the simplest way. Although the temperatures 
of the heater, cooler and regenerator vary with time because of 
their transient behaviour, they are assumed uniform along the 
different chambers. Only the expansion and compression space 
volume will vary with time and angle. The heater, cooler and 
regenerator volumes together with the swept and dead volumes 
will be considered as constants. 

Although the Schmidt analysis does not include the internal 
hydraulic or thermal losses among others irreversibilities, it will be 
precise enough to explain the transient Whispergen performance 
during the variation of the source temperature. The principal 
equations utilised in the Schmidt analysis are the expressions from 


(1) to (6). 


V E = V * E (1 cos 0 ) + V DE + V H 

(i) 

V c = ^(1 - cos (0 - a)) + V DC + V K 

(2) 

i p-sin a 
a = tan - 

t COS Oi 

(3) 


B = \Jt 2 + 2 tv cos a + v 2 



2 mRTc 

V se (S-Bc os (0 - a)) 



In equations (1) and (2), the terms Vde + V h and Vdc + Vk 
represent the total dead volumes from the heater and the cooler, 
respectively. The volume occupied by the nitrogen inside of the 
regenerator was neglected [9] in comparison with the other 
volumes. 

Fig. 2 illustrates the evolution of the pressure and the volume 
along the engine when the source temperature decreases. A simple 
analysis from a SE resolves that an increase in the source temper¬ 
ature represents an increment in the power output for a fixed sink 
temperature. The PV diagram for the three source temperatures 
illustrates in this figure that the work in the compression or cooler 
piston is maintained constant, while the work in the expansion or 
the heater piston is continuously reduced as the source tempera¬ 
ture decreases. 

By representing the different source temperatures or the 
exhaust gas temperatures against the experimental and theoretical 
power output, a lineal relation between both of them is found. This 
relation is true for both the Schmidt analysis and the Whispergen 
unit as Fig. 3 illustrates. Expression (7) can be used to calculate the 
power output from a Whispergen unit using the power obtained 
with the Schmidt's analysis at different heat source temperatures. 
This expression was experimentally validated in the Whispergen 
unit between 640 K and 720 K. 


^Whispergen [W] = Pgchmidt [W] - 3.54 -T H [I<] + 661.3 (7) 

Although this experimental expression could be used to calcu¬ 
late the power output, an expression, which reflects the influence 
of both the cold and hot source temperatures, is needed. 

The Beale number and the West number are constants that are 
experimentally determined 11 . They can be utilised to approxi¬ 
mately calculate the power output from a SE. The Beale number (8) 
estimates that the power output from a SE is proportional to the 
mean pressure, swept volume and speed. 



Pi 


Pm • ^se * N 



S 


t + 2tXp)£ + 


4 tX R 
1 + t 


+ v 4 s 2 Xdc 





s « « « 


Fig. 1. Stirling engine distribution and of four-cycle double acting alpha configuration. 
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Fig. 3. Comparison between Schmidt analysis and experimental data from Whisper- 
gen unit. 


power output from the engine by varying the hot source temper¬ 
ature in the combustion chamber. 

As explained previously, the power output from this Stirling 
engine is controlled via the airflow entering into the burner 
because the fuel pump introduces diesel into the combustion 
chamber for a control system. This control measures the oxygen 
content in the exhaust gases, and it keeps a prefix fuel-air ratio. 
Depending on this airflow, different hot source temperatures are 
reached. 

By varying the excitation power into the brushless alternator 
[29], the revolutions from the SE are practically kept constant 
during the running mode. This process provokes a constant thermal 
resistance and consequently, the heat flow, which is moved by the 
nitrogen inside the pistons between the hot and cold source, de¬ 
pends principally on the difference in the temperatures between 
the nitrogen in the cold side and hot side and the temperatures 
from the sink source and heat source. 

In the expansion space (hot side), the heat flow can be obtained 
with the expression (11): 


When William Beale performed the tests to determine this 
experimental number, he used similar source temperatures. For 
this reason, West overcame plus experiments where the heat 
source and the cold source temperature were varied. A similar 
relation was found, even though the heat and cold sources tem¬ 
peratures were also in the equation (9). 

w _ A _ R fT e + T c \ /qi 

S Pm-Vse-N-(T e -T c /T e + T c ) n \T e - Tq) y > 

For this experimental number, this dependence on the tem¬ 
perature term is implemented in the lump model, which improves 
the proposed model from the previous work of the authors. The 
output power is obtained with the expression (10) instead of using 
a straight line for the stopping process similar to the preceding 
model. This equation uses the temperature of the two lumps with 
two correction constants K\ and K 2 , which are experimentally 
obtained. 


<7 in = tt-ATh = T ( r H - Tu) (11) 

k H 

where Rh represents the thermal resistance from the hot side. 
Although this resistance depends on different factors (such as 
frequency, properties from the nitrogen, running temperatures or 
the heating method), its value can be considered constant [12] to 
simplify the SE analysis. In Ref. [9], this heat flow is calculated 
using the expression (12), where hu represents the convection 
coefficient. 


Qin = h H A„ (t h - Tu) (12) 

Following the same argument for the cold side, the heat flow 
that the nitrogen releases into the sink source can be calculated 
with the expression (13) or similarly (14 . The principal heat flows 
are depicted in Fig. 4. 


(,0) 

6. Lump model 

In a previous work from the authors, a transient model based on 
two lumps was developed within the TRNSYS simulation environ¬ 
ment [24]. In that article, the dynamic response from the cogene¬ 
ration model was principally reduced to the behaviour of two 
masses and one heat exchanger. The mass of these two lumps and 
the heat transfer coefficients were sized to obtain the temperature 
profiles of the exhaust gases and water outputs, respectively, when 
the sink temperature was varied. 

New experiments were carried out on the test bench, where the 
air mass flow entering into the combustion chamber was varied. A 
reduction in the air mass flow provokes the decline of the fuel 
entering into the burner and consequently, the heat source tem¬ 
perature is reduced during these tests. Comparing these new 
experimental results with the results obtained in the original lump 
model, important discrepancies were found, and an improved 
model was required. 

These new experiments were conducted by looking for a 
mechanism to vary the power output from the engine. This model 
of the cogeneration device allows us to indirectly control this 



Fig. 4. Temperatures and heat flows into the SE unit. 































A. Cacabelos et al. / Applied Thermal Engineering 73 (2014) 606-619 


611 



<jout = ^-AT C = ±(t l - T c ) (13) 


Qout = Me (Tl-Tc) (14) 

Consequently, the values of q in and q out will depend principally 
on the difference in the temperatures between the hot side tem¬ 
perature and the nitrogen temperature into the expansion space 
and on the difference in the temperatures between the sink tem¬ 
perature and the temperature from the nitrogen within the 
compression space as shown in Fig. 4. Different tests were per¬ 
formed in the cogeneration unit, and they showed that the engine 
operates at a similar speed during the starting and the steady state 
periods. However, during the stopping period, the revolutions from 
the engine suddenly drop. 

Therefore, the performance of the global system can be reduced 
to the evolution around the time of the two lumps as represented in 
Fig. 5. These lumps represent a simplified version of the SE, and 
they will replicate the behaviour of the real system by means of a 
method, which could be applied to other machines. 


The reduced equation (15) was used to represent the heat flow, 
which circulates between the lumps, and it uses an equivalent heat 
transfer coefficient. 


Qout — Q-lumps — (')lumps ^L2J 0^) 

Therefore, because the SE Whispergen unit runs with a similar 
rotating speed during the starting process and during the steady 
state period, it used the same heat transfer coefficient for these 
periods. By contrast, during the stopping process, it used a different 
coefficient of transference to take into account the sudden drop of 
the working speed and consequently the variation of the heat flow 
between the heat and cold sources. 

Although the speed reduction decreases the times per second 
that the nitrogen release or the capture heat from the sink or the 
heat source, respectively, the period of time to exchange this heat 
increases and provokes a bigger coefficient of transference. There¬ 
fore, the engine stops quickly when the fuel entering into the 
combustion chamber is cut off. 

The expression (16) is introduced in the TRNSYS simulation 
environment, and it controls the heat flow between the two lumps' 
model, which represent the heat exchange between the expansion 
chamber and the compression chamber from the Whispergen unit. 


Qi 


umps 


(Tu ~ T L 2l*<UA, lumpSi { ^<> 3 ° 

(Tu ~ T L2)*(UA) lumpS2 1 543 

0 in other case 


(16) 


This expression also contains information about the perfor¬ 
mance and control logic of the SE. The heat flow between the two 
lumps is allowed when different conditions are fulfilled. However, 
when the exhaust gas temperatures are greater than 140 °C, the 
engine starts to turn, and the heat exchange between the lumps is 
allowed. When there is no fuel going into the combustion chamber 
but the exhaust gases temperature is still over 270 °C, the heat 
exchange between both lumps is also allowed. These conditions 
correspond to the logical operation of the SE and are fulfilled in all 
of the undertaken experiments in the cogeneration unit. 

Fig. 6 illustrates the cogeneration unit that is simulated in the 
TRNSYS simulation environment. LI and L2 represent the lump 
units that simulate the behaviour of the SE, and ST_HE represents 
the heat exchanger placed in the exhaust gases. 



► 


Fig. 6. Schematic of the propose model in the TRNSYS simulation environment. 
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7. Heat exchanger arrangement into the cogeneration plant 

The original Whispergen unit uses an external heat exchanger to 
recover the heat flow within the exhaust gases. When the SE rea¬ 
ches the steady-state regime, the exhaust gases reach temperatures 
of approximately 400 °C. This heat flow represents an important 
quantity of energy compared with the energy introduced by means 
of the fuel inlet. For this reason, the Whispergen unit recovers this 
heat flow by means of a shell and tube heat exchanger using the 
cooling water (sink source) going out of the block and conse¬ 
quently, the exhaust gas temperatures decrease by approximately 
80 °C. 

Fig. 7 illustrates the original distribution of the heat exchanger 
within the cogeneration unit. A SE is primarily formed by three 
heat exchangers (heater, cooler and regenerator) that are care¬ 
fully designed to fit the compromise between different factors, 
such as pressure drop or heat exchange. The heat exchanger that 
corresponds to the regenerator was not represented in Fig. 7 
because this regeneration is considered perfect, and it corre¬ 
sponds with the internal performance of the SE that is not taken 
into account in this analysis. The equations (17)—(22) represent 
the SE behaviour and the principal exchanger energy inside the 
cogeneration system. 

QhE 1 — tHgas -Cp gas • (T L i — Tg as ) = n7 wa ter' water • (Tout — T L2 ) 

(17) 

nigas = m air + m fuel (18) 

JVfCn^p = fnjuei LHV - rh g c Pgas (T og - T lg ) 

-(Slumps(Ti ~ Pi ) - Power (19) 

MQ2-JP = ( UA) lumps (T.1 - Tn) - ™ water 1c p waK r(T.2 ~ Tret) 

( 20 ) 

Qcomb ~ ^fuel" LHV (21) 


Power = W ne t • n (22) 

The new configuration represented in Fig. 8 was simulated for 
the purpose of investigating the ETTR variation. This simulation 
was carried out in the TRNSYS simulation environment, and the 
obtained results were compared with the original data from the 
cogeneration plant that is represented in Fig. 7. By preheating the 
mass flow air that it is entering into the combustion burner, an 
important reduction in the fuel consumption can be obtained for 
the same power output. This reduction in fuel consumption is 
especially interesting when the system is tracking the electrical 
load. 

Equation (23) represents the equilibrium in the newly proposed 
heat exchanger. This final exchanger preheats the ambient air up to 
temperatures near the exhaust gases from the combustion cham¬ 
ber, and the UA from the heat exchanger was selected to maintain 
the same gas outlet temperature. 

Qhe 2 = thgas • Cp gas • (T L i — Lgas) = m air 'Cp air • (T in - T amb ) (23) 


8. Calibration methodology 

The capacitances from the block and from the burner, the 
coefficient of transference between the lumps, the coefficient of 
transference in the final heat exchanger and the constants k\ 
and k 2 from the power output equation were introduced as 
input variables in the model. They were obtained through an 
iterative process using TRNSYS and GenOpt. This optimisation 
program tries to minimise the objective cost function using a 
predefined algorithm. The selected cost function was the RMSE 
of the principal values from the SE: the power output, the 
exhaust gas temperatures from the burner and the temperature 
of the block. Because the input variables are continuous and the 
cost function is differentiable, the Hooke-Jeeves algorithm can 
be used to search for the best result. During the iterative pro¬ 
cess, all the input values were varied between acceptable 
values. 

This calibration process was performed using data from the tests 
with ninety percent of the air mass flow, and the model was 
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Fig. 7. Original heat exchanger distribution. 


Fig. 8. New proposal for the heat exchanger distribution inside the system. 
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validated with the experimental data from the tests with seventy, 
eighty and one hundred percent of the air mass flow. 

9. Results and discussion 

9.2. Experimental results 

For all of these dependencies and with the aim of validating the 
newly proposed model, different experiments were performed to 
collect the performance of the Stirling engine at different heat 
source temperatures. Because the sink source temperature depends 
principally on the heating system, all experiments were performed 
with a constant sink source temperature of 60 °C. 

The collected data from these experiments are represented in 
Figs. 9-11. In Fig. 9, it can be observed that the heat source tem¬ 
peratures in the steady state period for the different working points 
are in the range between 350 and 450 [°C]. 

During the first fifteen minutes, the Whispergen engine pre¬ 
sents the same behaviour and is independent of the set point for 
the airflow. During the starting process, the engine checks different 
sensors and parameters to ensure a good starting process, and it 
follows the same starting ramp until minute 14 is reached. 
Furthermore, once the exhaust gases from the burner achieve a 
temperature of 150 °C, the engine is forced to turn. 

For the power output, these differences are greater as shown in 
Fig. 10. At maximum power, the airflow is set at 100%, and the 
power output reaches values of approximately 1 kW. When the 
airflow is set to 70%, the power output in the steady state period has 
a value near 700 W. 

However, when the cogeneration system is forced to work with 
a 60% quantity of air, the power output is not sufficient to supply 
the consumption of energy from the test bench, and the engine has 
to stop to avoid the discharge of the batteries. For this reason, the 
collected data from the test with 60% air is represented, but it is not 
considered within the simulations. 

The voltage tracking from the system keeps the frequency from 
the alternator constant and consequently the speed from the en¬ 
gine. The brushless alternator from the Whispergen unit produces a 
three-phase alternating energy. The frequency of this current was 
measured by an oscilloscope during the experiments, and it was 
found that the values were in the range 150-170 Hz in the starting 
and steady-state periods. These values correspond to a speed be¬ 
tween 1500 and 1700 rpm for expression (24), where P represents 
the number of poles per phase. 
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Fig. 10. Power output working at different air mass flow. 
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P 


(24) 


The three-phase energy is converted into direct current with a 
rectifier, which assures the correct charging process for the battery 
bank at 24 V DC. Because the voltage from the system is maintained 
constant, the intensity going out of the SE follows the same varia¬ 
tion as the power output as shown in Fig. 11. 

Despite the differences for the source temperatures and the 
power output, the power efficiency from the engine approaches 
10%, and there is no clear disparity between the variety of running 
points. Due to the thermal inertia from the system, this efficiency is 
greater in the test of 60% and 70% air during the starting period, but 



0 1000 2000 3000 4000 5000 


Time [s] 


-•- Intensity with 60 % of air 

o Intensity with 70 % of air 

-▼-Intensity with 80 % of air 

-a - Intensity with 90 % of air 

— —■- — Intensity with 100 % of air 


Fig. 9. Exhaust gases temperatures running at different mass flow inlet. 


Fig. 11. Intensity from the alternator working at different source temperatures. 



























614 


A. Cacabelos et al. / Applied Thermal Engineering 73 (2014) 606—619 


>> 

o 

c 

0) 

o 


LU 



0 1000 2000 3000 4000 5000 


Time [s] 


-•- Electrical efficiency with 60 % of air mass flow 

o Electrical efficiency with 70 % of air mass flow 


-▼-Electrical efficiency with 80 % of air mass flow 

-a - Electrical efficiency with 90 % of air mass flow 

-■- — Electrical efficiency with 100 % of air mass flow 


Fig. 12. Electrical efficiencies at different air mass flows. 
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Fig. 14. Total efficiencies at different air mass flows. 


once the steady state is reached, these differences are reduced as 
Fig. 12 illustrates. As the air mass flow increases, the power effi¬ 
ciencies are only slightly larger, but these differences are not suf¬ 
ficient to extract a valid conclusion, and they cannot be used to vary 
the ETTR ratio. 

In Figs. 12 and 13, the obtained electrical and thermal effi¬ 
ciencies measured at different air mass flow inputs are illustrated. 
These experimental data show that the engine tends to work at a 
fixed working point that is independent of the source temperature. 
As the air mass flow decreases, the thermal efficiency slightly in¬ 
creases, but the power output drops below the critical value of 
600 W, where the SE is forced to stop. 

Figs. 13 and 14 show that after 30 min of running time, the 
cogeneration system starts to produce a considerable quantity of 


heat with a good efficiency. This transient performance shows the 
importance of developing transient models that are easy to inte¬ 
grate into simulation tools because the starting time is an impor¬ 
tant period compared with the total running time especially if the 
cogeneration system is not oversized, and it is turned on several 
times per day. 

For the stopping time, Fig. 15 and Table 2 represent the total 
energy produced during this downtime. This energy is also 
important to consider, especially during the first ten minutes, when 
the engine is still turning. 

The ETTR ratio from the simulated Whispergen unit represents 
the relation between the total power output and the total thermal 
energy obtained through the SE cogeneration device, and it can be 
calculated with expression (25). 
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Fig. 13. Thermal efficiencies at different air mass flows. 


Fig. 15. Electrical and thermal power for the different test during the stopping time. 
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Table 2 

Thermal and electrical power during the stopping time. 



60% Air 

70% Air 

80% Air 

90% Air 

100% Air 

Power [kWh] 

0.041 

0.052 

0.064 

0.068 

0.121 

Qth [kWh] 

3.28 

3.37 

3.17 

3.54 

4.75 


ETTR = p- (25) 

' th 

The total thermal energy P t h is measured by means of the 
circulating water circuit using the specific heat from the water, the 
mass flow and the difference in the temperatures between the inlet 
and the outlet as expressed in equation (26). 


P t h — tfiwater' c p wa tei ' (^out — ^ret) (26) 

Using both expressions, the ETTR ratio is calculated for all the 
different experiments, and it is illustrated in Fig. 16. This ratio is 
approximately 0.15 using 80-100% of the maximum mass flow air. 
As the mass flow is reduced, this ratio slightly decreases, even 
though no important variations are obtained. 

9.2. Calibration results 

Table 3 represents the obtained results in the optimisation 
process, where different inputs were varied for the calculated cost 
function considered in the calibration process. 

Fig. 17 represents the experimental and simulated data from 
the test with ninety percent of the air. In this figure, the output 
power and the exhaust gas temperatures at the burner output are 
depicted. In Fig. 18, the experimental and simulated data from the 
test with 70 percent of the air are also represented. This simulation 
was performed with the values of the capacitances and the co¬ 
efficients of transference that were obtained with the experi¬ 
mental data from the test with 90 percent of the air, and the 
obtained results show that the model adopted is valid with 
different air mass flows. 

Similar agreements are found within the tests of 80 and 100% 
of the air. Table 4 summarises the obtained RMSE from the 
different tests. It can be checked than the water outlet tempera¬ 
ture, the gas outlet temperature and the power output, which 
represent the most important outputs, are predicted with an er¬ 
ror less than 5%. 



Fig. 16. ETTR along the different test during steady state period. 


Table 3 

Optimisation results. 


Calibration results 


Capacitance_block 

54.24 

[kJ/K] 

Capacitance_burner 

12.86 

[1<J/K] 

Ualumps_starting-steady-state 

10.59 

[W/K] 

Ualumps_stopping 

15.65 

[W/I<] 

UA_final 

20.04 

[W/K] 

k\ 

0.19 

hi 

I<2 

20,838.69 

hi 
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Fig. 17. Experimental and simulated data from the power output and exhaust gases 
temperature at 90% air. 
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Fig. 18. Experimental and simulated data from the power output and exhaust gases 
temperature at 70% air. 



Table 4 

RMSE of the principal outputs from the different tests. 


% Air 

Gases 

Water 

Power 

70 

3.78% 

1.73% 

0.28% 

80 

3.95% 

2.32% 

0.37% 

90 

3.63% 

1.64% 

0.27% 

100 

3.88% 

1.78% 

0.30% 
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ETTR 80 



ETTR 90 



ETTR 100 




These results demonstrate that the values of the capacitances 
and the coefficients of transference that are obtained through the 
optimisation process are valid to simulate the behaviour of all 
possible transient behaviour from our SE. 

Fig. 19 represents the experimental and simulated ETTR ratios 
for the experiments with 80%, 90% and 100% of the air mass flow. 
These figures show a good agreement between the real and 
simulated ETTR data during the steady-state period. They are all 
approximately 0.15 as explained previously. 

9.3. Results of the modified SE cogeneration unit 

Fig. 20 represents the ETTR ratio from the original cogeneration 
plant, and it is compared with the ETTR ratio obtained in the 
simulation of the new design. This comparison was made for the 
test with 70% of the air mass flow. 

Similar results were obtained for the rest of the tests for the 
different air mass flows. These results are represented in Table 5. In 
the original distribution, the ETTR ratio is near 0.15 for all the ex¬ 
periments. However, with the new distribution, an ETTR ratio near 
0.26 is obtained for all the simulated cases. The simulations were 
performed with the same mass fuel inputs obtained from all the 
experiments and that correspond to each air mass flow rate. 

Figs. 21-24 represent the obtained results for the different 
analysed cases. In these figures, the simulated power of the original 
distribution is compared with the simulated power from the new 
arrangement. Moreover, the principal outlet temperatures are also 
compared. 

The values corresponding to the steady-state values for the 
original distributions and the new distributions are summarised in 

Table 6. 

In general, due to the increase of temperature in the exhaust 
gases from the combustion chamber, the temperature in the 
expansion space also increases and corresponds with an increment 
for the difference in the temperatures between the cold source and 


Table 5 

ETTR from the original and the new heat exchanger distribution. 



% Air 

ETTR_original 

ETTR_new 

Original distribution 

70 

0.143 

0.257 


80 

0.147 

0.265 


90 

0.148 

0.268 


100 

0.148 

0.264 


Fig. 19. Experimental and simulated ETTR from test with 80%, 90% and 100% of air. 
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Exhaust gases temperature and power output with 70 % of air 



Simulated power of the original SE 
a Simulated power of the SE with the new HE 
Simulated outlet temperature of the original SE 
o Simulated outlet temperature of the SE with the new HE 


Exhaust gases and water ouput temperatures with 70 % of air 



Time [h] 


—•- Exhaust gases temperature of the original SE 

o Exhaust gases temperature of the SE with the new HE 

—▼-Output water temperature of the original SE 

- Output water temperature of the SE with the new HE 


Fig. 21. Power output and principal outlet temperatures with 70% of air mass flow rate. 


the heat source, which directly implies an augmentation on the 
theoretical efficiency corresponding to Carnot's theorem. 

The security system from the Whispergen unit does not allow an 
increase in the exhaust gases temperature of over 480 °C, and the 
engine is forced to stop if this happens. This protocol assures the 
preservation of the sealing and assures that the working gas is kept 
inside the cylinders. 

Table 6 shows that with the new distribution, a power of 
approximately 1 kW can be reached with 70% of the air mass flow 
and supposes a 30% reduction in the fuel consumption to produce 
the same power output as 100% of the air mass flow in the original 
SE distribution. 

These results reveal that an important fuel saving could be ob¬ 
tained with the new distribution in the cogeneration plant. This 
new arrangement becomes especially interesting if the cogenera¬ 
tion unit is for an electrical load or can toggle between the heating 
and charging mode similar to the Whispergen unit. 


Exhaust gases temperature and power output with 80 % of air 



-±— Simulated power of the original SE 
a Simulated power of the SE with the new HE 
Simulated outlet temperature of the original SE 
o Simulated outlet temperature of the SE with the new HE 


It is important to remember that this is a preliminary analysis, 
where the different pressure losses or the new heat exchanger 
parameters corresponding to the nitrogen working at higher tem¬ 
peratures were not taken into account. 

Other heat exchanger distributions similar to Fig. 25 can obtain 
an intermediate ETTR ratio if desired. Furthermore, this ETTR ratio 
could easily be made variable with an adequate system of linkage, 
which could track the different loads that increase the dynamic 
response from the model. 

10. Conclusions 

In this article, the performance of a cogeneration unit that is 
principally composed of a SE, an alternator and a heat exchanger is 
studied. The obtained experimental results that work at different 
mass fuel inputs are presented and analysed. These tests corre¬ 
spond to a variation in the source temperature. 

Exhaust gases and water ouput temperatures with 80 % of air 



—•- Exhaust gases temperature of the original SE 

o Exhaust gases temperature of the SE with the new HE 

-Output water temperature of the original SE 

_a- Output water temperature of the SE with the new HE 


Fig. 22. Power output and principal outlet temperatures with 80% of air mass flow rate. 
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Exhaust gases temperature and power output with 90 % of air 


Exhaust gases and water ouput temperatures with 90 % of air 
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-r-Output water temperature of the original SE 

-a - Output water temperature of the SE with the new HE 


Fig. 23. Power output and principal outlet temperatures with 90% of air mass flow rate. 


Exhaust gases temperature and power output with 100 % of air 
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Fig. 24. Power output and principal outlet temperatures with 100% of air mass flow rate. 


A new transient model of a Stirling engine based on two lumps 
and one heat exchanger is presented. This model is adapted to the 
different transient states that correspond to different air mass flow 
inputs. 


Table 6 

Steady state values for the most important variables. 



% Air 

m fU ei [kg/h] 

Tom ra 

Texh ra 

Tgas [°C] 

Power [W] 

Original 

70 

0.56 

64.73 

72.49 

364.60 

691.33 

distribution 

80 

0.64 

68.35 

81.45 

397.96 

795.33 


90 

0.73 

69.97 

89.03 

420.16 

876.95 


100 

0.81 

71.26 

96.81 

455.04 

993.43 

New 

70 

0.56 

62.30 

67.71 

464.27 

1030.11 

distribution 

80 

0.64 

65.48 

79.97 

520.46 

1177.42 


90 

0.73 

66.21 

87.89 

555.34 

1274.58 


100 

0.81 

67.48 

100.38 

623.20 

1443.44 


Furthermore, the starting and stopping periods are better 
simulated than the new relations founded in the model. The heat 
flow between the different lumps was also accommodated to the 
performance of the model, and the different heat transfer co¬ 
efficients were calculated for the engine speed. 

The data obtained from the experiments were used to calibrate 
the model and had a RMSE value below 0.4% for the generated 
power. 

The new model has been used to investigate the behaviour of 
the SE when the thermal energy from the exhaust gases is used to 
preheat the air in the combustion chamber. This mechanism allows 
the increase of the electrical efficiency and an important fuel saving 
when the system tracks the power consumption. 

The air mass flow rate entering into the burner was also varied, 
and the real power output from the engine was analysed together 
with the electrical and thermal efficiencies. 
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Fig. 25. Intermediate solution for the heat exchanger distribution inside the cogeneration plant. 


Additionally, the EUR ratio was highly investigated in this 
commercial SE cogeneration device. Important variations in this 
EUR ratio were obtained by modifying the heat exchanger 
arrangement along the cogeneration plant. The new design can 
reach an EUR ratio of 0.26 in comparison to 0.15 from the original 
scheme. 

For the produced electricity, 30% of fuel saving can be achieved. 
Moreover, if the unit it is adapted to work at higher pressures and 
temperatures, power over 1.5 kW can be obtained in the alternator. 
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